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INTRODUCTION 


The  use  of  minimum  mean-square  error  (MMSE)  filtering  for  prediction,  estimation 
and  smoothing  has  grown  in  popularity  during  the  last  few  years.  To  a large  degree  this 
interest  may  he  attributed  to  the  corresponding  growth  in  the  areas  of  real-time  computing 
hardware  and  software.  The  degree  of  sophistication  reached  in  these  areas  has  made  it  feasi- 
ble to  construct  dedicated  hardware  componentry  or  interactive  software  algorithms  capable 
of  implementing  MMSh  filtering. 

As  originally  developed  by  Wiener,*  the  theoretical  techniques  of  MMSH  filtering  per- 
mitted us  to  solve  analytically  for  the  linear  filter  impulse  response.  This  would  minimize 
the  mean-squared  error  between  the  filter’s  output  signal  and  a desired  signal.  We  have  con- 
siderable freedom  in  choosing  this  desired  signal:  the  specific  application  determines  its  choice. 

One  possibility  enables  us  to  use  the  output  of  a linear  filter  and  to  predict  the  future 
values  of  a time  series  when  we  know  its  present  and  past  values.  The  optimal  filter  impulse 
response  (which  minimizes  the  mean  square  of  the  error  between  the  filter  output  and  the 
selected  time  series  value)  appears  (in  the  continuous  case)  as  the  solution  to  a Fredholm 
integral  equation  of  the  first  kind.  Solution  strategies  for  this  equation  have  been  developed 

by  a number  of  authors ‘+  for  several  signal  and  noise  situations.  For  the  case  of  discrete 

data,  the  solution  for  the  optimal  L-length  filter  coefficients  appears  as  the  solution  vector, 
w ■*.  to  the  well-known  discrete  Wiener-Hopf  matrix  equation: 

Rh=PA.  <11 

In  Equation  1 , R is  the  L X L.  autocorrelation  matrix  of  the  input  time  series,  h is  the  L X 1 
vector  of  the  optimal  weight  coefficients  (the  discrete  impulse  response  of  the  filter)  and  PA 
is  the  L X I autocorrelation  vector  with  autocorrelation  elements 


^xx* A) 

0XX(A  + *> 


XX*  A + L " 1 ) J 

where  A is  the  prediction  distance  in  sampling  intervals. 

The  main  difficulty  in  solving  Equation  1 is  the  estimation  of  the  true  autocorrelation 
lags.  0xx(k ).  lor  the  time  series  under  examination.  This  in  itself  is  an  area  of  active  current 
research.  In  actual  systems,  we  are  given  a finite  amount  of  data  from  which  to  estimate  the 
actual  lag  values.  This  introduces  some  degree  of  uncertainty.  This  work  examines  the  theo- 
retical case  of  assuming  the  autocorrelation  lags  known  exactly  and  develops  an  analytical 
solution  to  Equation  1 . Solutions  developed  for  this  ideal  case  are  useful  for  comparison  with 
actual  processors  which  approximate  solutions  to  h from  finite  segments  of  noisy  data. 
Although  this  is  a natural  application  of  the  analysis  in  this  paper,  we  will  confine  the  present 
work  to  examining  the  properties  of  the  theoretical  solutions  to  Equation  I . 

The  specific  case  to  be  examined  is  that  of  bandlimited  (BE)  signals  corrupted  by 
lowpass  noise.  An  example  is  the  signal  generation  by  a Gaussian  random  process  of  a band- 
limited  signal  and  the  noise  characteristics  of  the  media  exhibit  the  familiar  “IF"  frequency 


t 


j 
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absorption.  An  analytical  technique  known  as  the  method  of' undetermined  coefficients  is 
applied  toward  the  solution  of  Equation  1 for  this  case.  Early  theoretical  work  done  by 
Zadeh  and  Ragazzini-  devised  this  method  for  continuous  systems.  The  method  was 
adapted  to  the  discrete  case  by  Krut’ko''  and  Solodovnikov.9  More  recently,  general 
properties  of  the  discrete  method  were  examined  by  Satorius  and  Zeidler.^  The  discrete 
method  was  then  applied  by  Zeidler  et  al°  to  solving  Equation  1 for  the  case  of  multiple 
sinusoids  in  uncorrelated  noise  and  by  Satorius  and  Zeidler9  to  the  case  of  multiple  sinu- 
soids in  lowpass  noise.  To  fully  examine  the  effects  of  a nonzero  bandwidth  signal  upon  the 
resulting  optimal  filter  structure,  the  input  is  limited  to  one  complex  bandlimited  signal  in 
lowpass  noise.  Motivation  for  examining  this  case  stems  from  Reference  8.  That  work 
reveals  that  as  long  as  the  frequencies  of  the  complex  exponentials  comprising  sinusoids 
are  sufficiently  separated,  there  is  effectively  no  interference  between  each  complex  expo- 
nential solution.  Solutions  corresponding  to  real  sinusoids  then  become  linear  superposi- 
tions of  the  complex  exponential  solutions. 


ANALYSIS 

The  matrix  Equation  1 may  be  expanded  into  its  L components  to  become 


L-l 

^ 0xx(k)  h(k)  = 0xx(t  + A):  C = 0,  1 L-l.  (3) 

k=0 


For  the  specific  case  of  one  BL  signal  in  lowpass  noise,  the  form  of  the  autocorrelation  lag 
0XX(C)  becomes 


0XX(£)  = °jvj  e^NlU  + a-  e^’sl^  e-^s® 


(4) 


In  Equation  4, 

i 

o“  = noise  mean-square  power 
= signal  mean-square  power 

as  = correlation  parameters  of  noise  and  signal,  respectively 
cos  = radian  center  frequency  (relative  to  Nyquist)  of  complex  BL  signal 


Following  the  method  outlined  in  Reference  7.  we  assume  the  general  solution  to 
Equation  3 for  the  case  of  a complex  BL  signal  in  lowpass  noise  to  be: 

h(k)=  Bm  z^  + C,  6(k)  + C2  5(k  - L-  1),  k = 0.  1 L-l.  (5) 

m=  I 

In  Equation  5.  the  Bm.  C)  and  C'2  are  complex  coefficients,  as  yet  undetermined:  the  zm  are 
damped  complex  exponentials  of  the  form  zm  = exp(pm  +jdm):  and  the  delta  functions 
6(  • ) are  due  to  the  end  effects  of  the  finite  length  filter. 

The  solution  strategy  is  to  substitute  Equation  5 for  the  optimal  weight  vector  and 
Equation  4 for  the  autocorrelation  function  into  Equation  3 and  solve  the  set  of  resulting 
equations  for  the  unknowns  Bm,  zm.  (']  and  C'2-  Applying  this  technique  yields  six 
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equations  in  six  unknowns,  which  may  then  he  partitioned  according  to  the  following  con- 
sistency requirements: 

1.  The  coefficients  of  in  the  resultant  substitution  must  he  equal  for  every  f of 
the  equation; 

2.  I’lie  coefficients  of  e-0^  and  eQN^  must  similarly  equate  and 

3.  The  coefficients  of  e,-tts  + and  e,us  + i^s1^  must  also  equate 


Condition  1 leads  to  a complex  quadratic  equation  which  may  he  solved  for  the  com- 
plex exponentials  /)  and  z\  Conditions  2 and  3 give  a set  of  four  equations  in  the  remaining 
four  unknowns  < B ] . B^.  I | , Cq.  Thus  the  unknown  parameters  of  Equation  5 may  he 
found  and  the  resulting  optimal  weight  vector  determined. 

Applying  condition  I gives  the  following  equation: 


N 


1 _ y C-«N 
I 'lU  1 


+ o; 


Z e-os-JWs  | _7  e«s-JWs 

m ' in 


= 0 


which  after  a few  operations  leads  to 

1 T ■ "| 

sinh  - 2 cosh  as  eJUJs  /m  + e|-t°s) 


+ a~  eJCJs  sinh  as(zj^  - 2 cosh  zm  + I ) = 0. 


(6) 


(7) 


Equation  A17  in  the  Appendix  shows  that  this  is  equivalent  to  the  condition 
N(zm)  = 0 (8) 


where  Ntz)  is  the  /-domain  polynomial  representing  the  numerator  of  the  power  spectral 
density  Sxx(z)  given  hy 


c . Ntz) _ 
•xx* n D(z)  1 


{0XX(£)} 


(d) 


In  Equation  L>.  I)(z)  represents  the  denominator  z-polynomial,  </>xx(C)  is  the  autocorrelation 
function  from  Equation  4 and  Z{* } signifies  the  /.-transform  operation.  Equation  8 states 
that  the  zm  are  located  at  the  zeroes  of  the  input  spectral  density . For  the  case  of  one  BL 
complex  line  the  zeroes  are  given  by  the  pair  {/| . Z2}  which  are  reciprocal  to  the  unit  circle 
on  tile  frequency  radial  0 | : 


z j = e_#i  e ^ 1 ( 10a) 

z.s  = e^e^l.  (10b) 


In  Equation  10.  /u  is  a measure  of  the  closeness  of  the  zero  to  the  unit  circle  and  for  the 
present  case,  /a  >«s- 

Applying  the  consistency  conditions  2 and  3,  we  obtain  the  following  set  of  equations: 
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= 0. 


The  substitutions  from  Equation  10  are  then  used  in  Equations  1 I from  which  the 
unknown  coefficients  may  be  determined  and  the  solution  for  the  weight  vector  becomes 

h(k)  = B,  e-^ke^lk  + B:e^kej^lk  + C,  5<k ) + Cs  5<k  - L + 1). 

k = 0,  1 L-  1.  ,12) 


ASYMPTOTIC  PROPERTIES 

Equation  12,  with  (B),  B2,C|.C2}  determined  from  Equations  1 I , gives  the  gen- 
eral solution  for  the  weight  vector  for  arbitrary  filter  length  L.  signal  bandwidth  parameter 
o fs  and  noise  bandwidth  parameter  ajs|.  However,  more  insight  into  the  properties  of  the 
optimal  impulse  response  may  be  gleaned  by  considering  a few  asymptotic  properties  of  the 
general  solution. 


LARGE  FILTER  LENGTH.  L 

First,  consider  the  form  of  the  impulse  response  as  the  filter  length  L becomes  larger 
and  larger  ( L — °°  in  the  limit).  Note  that  the  term  e/jk  present  in  Equation  1 2 is  unbounded 
for  increasing  k.  For  L unbounded,  this  will  cause  an  unbounded  solution  for  the  weight 
vector.  Thus,  if  we  consider  the  case  of  long  filter  length  L.  < L — °°  in  the  limit ).  the  coeffi- 
cients B2  must  be  equal  identically  to  zero  to  give  a stable  solution  for  htk).  With  this  sim- 
plifying assumption,  the  Equations  I la  and  1 Ic  become 


which  may  he  solved  for  B|  and  Cj  to  give: 

| | _ ea’N  - A*  ej0  1 | | i _ eas  - P e.i( » 1 - ws)  | 


B,  = 


C,  =e"aNA 


e«N  ~ A<  ej0  ! _ e«s  - V ej(0 1 " ws) 

| | _ e«s  "M  e.i(0 1 - ws>]  [e"aNA  _ e(-as  + jcos)A| 


e^NA_e(-as  + jws|Aj  ()4a) 


e“N  - A«  I - e“s  “ F c'-i'y  1 - wS) 

Making  the  approximation  B^  = 0 in  Ki] nation  I 1 b gives  the  following  expression  for  ('■>: 
e“  e)0 1 L 


(14b) 


C\  = e"“N 


(15) 


1 _ e-GN  - M 

But  for  the  long  filter  length  ease  ( L > 4//a),  the  factor  e~^k  asymptotically  approaches  zero 
and  thus  C"2  vanishes. 

Therefore,  for  the  case  of  very  long  filter  length  the  optimal  weight  vector  reduces  to 
the  following  approximation: 


h(k)  = B , e“^kej0lk  + C,  6tk) 


k = 0. 


L - 1 : L -*•  oo. 


( 16) 


The  weight  vector  given  by  Equation  1 6 has  two  parts:  a damped  complex  exponen- 
tial oscillating  at  the  frequency  of  the  spectral  zero  of  the  input  process,  and  an  impulse  of 
strength  Cj , existing  in  the  first  weight.  The  damped  exponential  decays  away  from  the 
beginning  of  the  filter  with  an  envelope  determined  by  the  value  of  jli.  the  zero-damping 
parameter.  For  very  small  ja,  the  oscillatory  function  is  almost  a pure  sinusoidal  term:  for 
larger  /a  the  solution  decays  to  zero  quickly.  The  impulse  C|6(k)  gives  a value  C'|  to  the  first 
(k  = 0)  weight  only.  Physically,  the  BL  signal  may  be  thought  of  as  causing  the  damped 
exponential  portion  of  hlk),  while  the  lowpass  noise  contributes  the  impulse  with  magnitude 
C|S(k). 

Once  the  impulse  response  hlk)  has  been  derived,  the  transfer  function  H(jco)  of  the 
optimal  filter  may  be  determined  by  taking  the  z-transform  of  h(kl  and  evaluating  around  the 
unit  circle  at  z.  = e-|aj.  Thus  from  Equation  I 6, 


L-l 


H(z)  = Z { h< k ) } = 


h(  k ) z 


-k 


= 2 

k=0 


L- , 

= B,  £ [e-^  + '0l  z-'|k+C, 
k=0 


or 


H(z) = B, 


_e-/aL  +j0\  L Z-L 


+ Ci  • 


(17) 


I _ 0-A<  + 1 z- 

Since  L -*■  °°.  the  exponential  e~^k  approaches  zero.  LIsing  this  approximation  and 
evaluating  at  z = e)co.  the  transfer  function  H(jco)  becomes 
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Thus,  the  optimal  transfer  function  has  a constant  background  component  of  complex  value 
C'l  (given  by  Equation  14b)  plus  a component  due  to  the  HE  signal  presence.  The  spectral 
peak  of  this  latter  term  is,  in  general,  a function  of  all  the  parameters  (as  seen  from  Equation 
14a)  and  is  not  examined  further  in  this  paper.  However,  an  enlightening  simplication  results 
from  considering  the  BL  signal  in  white  noise,  which  we  do  next. 
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An  important  property  of  the  impulse  response  tor  the  BL  signal  in  white  noise  is 
seen  from  liquation  24.  Namely,  the  oscillatory  portion  of  litkl  is  at  exactly  the  center  fre- 
quency. cos.  ot  the  signal  regardless  of  bandwidth.  Moreover,  the  impulse-type  noise  contri- 
bution due  to  ('  | has  disappeared  for  the  limiting  case  of  white  noise,  f urther  characteristics 
of  the  solution  are  seen  In  considering  the  effects  of  00  upon  B | . This  may  be  exam- 
ined by  first  multiplying  the  numerator  and  denominator  of  liquation  14a  by  e~QN,  giving: 

| C-«N  _ ejtf  1 | | | _ e«s  ~ A4  ej( 0 l - u>s)  | 

1 ej0  I _ C-«N  + «s  - A*  0j( 0 1 - ccs) 


|e-«NA  _0(-as  + jccs)A| 

C 


Taking  the  limit  as  cv^-  -*• 

fun  B,  e-^e-^l  |l  - e“s  " A4  LJ("  1 ’“si]  +.icos)A 

1 = . (26) 

aN  -*  00  e-*1  eJ°  1 

But  U i = tcs  for  the  white  noise  case  anil  thus 

B | = e-as^  ( | _ tf«s  - A4 , q.i^sA  ,27) 

B | is  thus  complex  with  a magnitude  determined  by  the  bandwidth  of  the  signal  (due  to  as) 
and  delay,  A.  and  with  an  initial  phase  determined  by  the  signal  center  frequency  cus  and  delay. 
When  liquation  24  is  substituted  into  liquation  27,  the  solution  for  the  weight  vector  becomes 

111  k 1 = t I - e°s  " F ) e_^k  e-'^s*  k + A ) (28 ) 


0 < k < L - I . 

The  solution  is  thus  seen  to  be  a damped  complex  exponential  with  initial  phase  tosA  and 
initial  magnitude  e-as^  ( 1 - eas  ~ For  increasing  A.  the  magnitude  decreases  as  e_as^: 
additionally,  the  initial  phase  increases  as  A is  increased. 

The  transfer  function  is  easily  obtained  from  liquation  IX  using  liquation  27  for  Bj 
and  the  fact  luat  C]  = 0 for  white  noise: 

0(-«s  +jcjs)A  ( , _ e«s  - A4, 

H(  icci  = : •.  (29) 

] _ e-A*  c.Kws  ~ w) 


SEPARATION  OF  SIGNALS  FROM  NOISE  USING  DELAY  PARAMETER.  A 

In  discrete  Wiener  prediction  theory,  we  estimate  the  value  of  a discrete  time  series  A 
samples  into  the  future  based  on  the  knowledge  of  its  present  and  past  values.  That  is.  based 
on  a knowledge  of  the  sequence*  { x( k ),  x( k — 1 ) x(k  - L + 1)},  we  form  an  MMSE  pre- 

diction ol  x(k  + At.  I h is  is  equivalent  also  to  predicting  the  present  value  of  the  time  series 

x( k 1 based  01  the  delayed  sequence  4 x(k  - A),  x(k  - A - I ) x< k - A - L + It}.  In  this 

section  we  show  that  we  are  able  to  separate  the  BL  signal  from  the  lowpass  noise  background 
using  a variable  prediction  distance  A.  A requirement  is  some  a priori  knowledge  of  the  ex- 
pected signal  and  noise  correlation  distances  (or.  equivalently,  signal  and  noise  3 cl B bandwidths). 

for  tillers  not  const  ruined  to  he  linite  length,  tins  sequence  extends  an  infinite  distance  in  the  past: 
x(  k I.  x(  k - 


Ui. 
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Using  liquations  14  uiui  1(>  we  see  that  the  expression  for  the  impulse  response  for 
long  filter  lengths  may  he  rewritten  to  retlect  the  dependence  of  the  B|  and  C]  coefficients 
upon  A: 

hA(k)  = B,(  A)  e“^k  eJy  lk  + C,(A)  6(k)  (30) 

k = 0 L - 1 : L 

where  B]l  A)  and  (']( A)  are  given  by  Equations  14a  and  14b.  respectively.  From  Equations 
14  we  see  that  for  fixed  signal  and  noise  parameters  (a\.  as,  0 | . cos.  p ) the  coefficients 
Bit  A)  and  C'|(  A)  (and  hence  the  filter  impulse  response)  are  functions  of  the  prediction  dis- 
tance A.  Since  B j ( A ) is  associated  with  the  BL  signal  characteristics  of  the  optima)  filter 
and  C'  j ( A)  is  associated  with  the  noise  contribution,  we  can  enhance  the  BL  signal  properties 
by  choosing  a value  of  A for  which  B|(  A)  dominates  Cj(  A). 

This  can  be  seen  clearly  in  the  following  example.  Consider  the  input  power  spec- 
trum as  shown  in  Figure  1 and  suppose  we  desire  to  use  the  MMSH  filter  impulse  response  to 
estimate  the  center  frequency  of  the  BL  signal.  For  the  parameters  given  in  Figure  I . the 
MMSH  filter  has  the  impulse  response  given  by 

hA(k)  = B,(A)  [.772  4t>0,7r|k  +C,(A)  6(kt  (31  I 


where  from  Equation  7.  the  / = / 1 . / s are  found  to  he 


z,  = .772  ej'-490,7r 


(32a) 


z2  = 1.295  e-K-490)7r.  (32b) 

From  Equation  31  we  can  see  that  the  impulse  response  now  becomes  a function  of  A via 
the  B |(  A)  and  C |(  A)  dependence.  Figure  2 shows  the  relative  magnitudes  of  B ]( A)  and 
C ((A)  as  a function  of  A for  the  signal  and  noise  parameters  given  in  Figure  1 . From  Fig- 
ure 2 we  see  that  for  small  values  of  delay  ( A < 7).  the  noise  coefficient  C']<  A)  is  much  larger 
in  magnitude  than  the  signal  coefficient  B](A).  The  resulting  impulse  response  would  then 
be  strongly  influenced  by  the  noise  component,  producing  a strong  impulse  for  the  k = 0 
weight.  This  is  shown  explicitly  in  Figure  3a.  which  is  the  real  part  of  the  complex  weight 
vector  for  the  value  A = I . There  is  a strong  impulsive  value  at  k = 0,  after  which  the  signal 
components  take  the  form  of  a damped  complex  exponential.  Note,  however,  that  for  A > 7 
in  Figure  2 the  B ] ( A)  coefficient  is  the  larger,  which  signifies  the  signal  components  should 
dominate  the  impulse  response.  This  is  clearly  seen  in  Figure  3b,  which  illustrates  the  almost 

1 complete  disappearance  of  the  impulse-type  value  at  k = 0 for  A = 7.  The  weight  vector  is 

very  nearly  a noise-free  damped  complex  exponential  which  should  lead  to  good  frequency 
estimation  properties. 

Figure  4 then  presents  the  magnitude  square  of  the  transfer  functions.  |H(u>)|-,  of 
the  weight  vectors  shown  in  Figure  3.  For  A = I (Figure  4a)  the  transfer  function  is  very 
much  dominated  by  the  noise  component  and  a highly  inaccurate  estimation  of  the  BL  signal 
results.  As  A is  increased  to  A = 7.  it  is  possible  to  estimate  frequency  better  from  the  trans- 
fer function  (Figure  4b).  One  effect,  however,  (seen  from  Figure  4b)  is  that  the  resultant  BL 

! signal  appearing  in  the  transfer  function  is  broader  than  the  original  BL  signal  shown  in  Fig- 

ure 1 . However,  the  spectral  peak  of  the  transfer  function  is  at  the  correct  center  frequency 
u)s  = ,5Cbr  of  the  BL  signal. 
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By  increasing  A to  a value  greater  than  the  correlation  distance  of  the  lowpass  noise, 
the  effects  of  the  noise  upon  the  MMSK  filter  response  may  be  diminished  considerably. 
However,  as  long  as  A is  not  increased  past  the  correlation  distance  of  the  Bl.  signal,  the  filter 
response  still  retains  good  characteristics  for  detecting  and  estimating  the  center  frequency 
of  the  BL  signal.  If  we  define  toy  as  the  3 d B bandwidth  of  a BL  process,  coy  = 2a  leads  to  a 
process  autocorrelation  function  of  the  form  0(  d)  = exp[-  a|C|  | . This  leads  to  defining  the 
correlation  distance,  I\  of  the  process  as 

r = 1/a,  (33) 

so  that  when  C = T,  the  value  of  </>(  (!)  = e“* . For  the  lowpass  noise  process  of  Figure  1 . the 
correlation  distance  of  the  noise,  F)nj,  then  is  given  by 

Fjyj  = 1 /cejsj  = 2.5  samples, 

whereas  the  correlation  distance  of  the  BL  signal,  Fs,  is  much  longer  and  given  by 
Ts  = l/as  = 20  samples. 

Thus,  the  value  of  A = 7 samples  chosen  from  Figure  2 falls  safely  within  the  range 
Fjyj  < A < Ts  required  to  give  good  signal  characteristics  in  the  transfer  function. 
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APPENDIX 

AUTOCORRELATION  FOR  CORRELATED  NOISE  AND  BL  COMPLEX  LINE 

The  autocorrelation  function  0(C)  and  power  spectral  density  Site)  are  related  by 
the  inversion  integral 


0(d ) = — I S(co)  dw  . 


Now  consider  the  discrete  form  S(z)  related  to  S(c o)  through  the  transformation 
(for  T = 1 ) / = ej°T  The  new  parameters  in  Equation  A1  may  be  found  by  the  following 
relations: 

dz  = j dec  -*  dco  = -j  /■'  dz 

co  = -7T  z = e-’^ 

co  = 7T  z = e-^  . 

The  last  relation  implies  the  path  of  integration  in  the  z-plane  is  counterclockwise 
around  the  unit-circle.  Thus,  in  the  z-plane  Equation  A1  becomes 


0(C)  = 37T  f s<z)  z* 


ResS(z)z  poles  of  S ( z ) 

Assume  first  that  S(z)  is  the  lowpass  noise  spectrum  given  by  S^jq(z): 


SNN,Z)  = 


- "Om  . . 

2 e 1 sinh 

(r'-e^T 


Applying  the  inversion  integral  (Equation  A2).  the  autocorrelation  function 
becomes : 


l2  ovi~  sinh  Oki  ) zA 

*NNltl  ~ (j*:, ) ‘ m 8 


Evaluating  Equation  A4: 

0NN,C)  = °n"  e N • c > 0 • 

A similar  evaluation  for  C < 0 gives 

0NN,C>  = °N2  °aN  C • 6<  0 • 

The  combination  of  Equations  A5  and  Ah  give 

1 -Ofv  I8i 

0NN<Ct  = aN~  e 


J 


r 


which  is  the  desired  autocorrelation  function. 

Next  consider  the  complex  BL  signal  with  the  spectrum  Sss(z)  given  by: 


Sss<z> 


1 2 ~“s  • 

2 as~  e 


sinh  a.. 


(z_c-“s  + )“s)  (z-l.c-“s -■!“»)  ' 

Again  applying  the  inversion  integral  gives  the  autocorrelation  function  0SS(8): 


(A8) 


(2  os~  sinh  as)  z®  e"*Ws 


, £>0 


( A9) 


Upon  evaluating.  Equation  A 9 reduces  to 


->  -a  8 jto,.  8 
«>ss  <»>  = Ve  S e s 

A similar  solution  for  8 < 0 gives 

->  a 8 jo:  8 

0SS(8)  = as-e  s e s . 


, 8>  0 . 


8 < 0 


(AIO) 


( A I I ) 


and  combining  Equations  AIO  and  A1 1 gives  the  required  autocorrelation  function  for  the 
BL  signal: 


0ss(8l  = os-  e 


-<xs  18 1 jtJs8 


(A12) 


The  linearity  of  autocorrelation  functions  now  allows  us  to  find  the  autocorrelation  t unc- 
tion. 0XX(8),  of  an  input  process  consisting  of  a complex  BL  signal  in  lowpass  noise  by 
simple  addition  of  Equations  A7  and  A 12: 

0XX(8)  = *ss(E>  + *NN(6)  ■ 

The  linearity  of  the  z-transform  operator  then  gives 


sxx(z)  ~ SNN(z)  + Sss(z)  ■ 
Performing  the  addition  in  Equation  A13  gives 
- _“N 


( A 1 3 ) 


2 a 


SXx*z> 


N 


sinh  a^t 


1 -«s  . , 
a - e sinh  cv 


After  considerable  algebra,  the  final  form  of  the  input  spectrum  becomes 
( - 2 z ) N</.) 


Sxx(z) 


l)(z) 


( A 1 4 1 


( A 1 5 1 
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who  re 


N(z)  = 

•>  . , / •>  iwt  j2to  \ 

°N~  smh  1 z~  - 2/  e s cosh  qs  + e s I 

1 iws  / i \ 

+ os“  e s sinli  crs  1 /-  -2/  cosh  + I 1 

(A16) 

IX  z.)  = 

(/-•'“*)  I,.,-*'!".)  (,.,1'K) 

( A 1 7) 

The  poles  of  Sxx(z)  are  the  roots  of  the  Equation  D (z)  = 0 and  are  the  pole-pairs  of  each 
separate  spectrum  from  Equations  A3  and  A8.  The  roots  of  the  Equation  N (z)  are  the  non 
tiivial  zeros  ol  Sxx(z).  Additionally,  there  is  a zero  at  the  origin  which  does  not  affect 
computation  of  tiie  spectrum  or  the  solution  for  the  weight  vector. 
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